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An inelastic neutron scattering study of the spin waves corresponding to the stripe antiferromag¬ 
netic order in insulating Rbo. 8 Fe 1 . 5 S 2 throughout the Brillouin zone is reported. The spin wave 
spectra are well described by a Heisenberg Hamiltonian with anisotropic in-plane exchange interac¬ 
tions. Integrating the ordered moment and the spin fluctuations results in a total moment squared of 
27.6 ±4.2/i^/Fe, consistent with S«2. Unlike XFe 2 As 2 (X = Ca, Sr, and Ba), where the itinerant 
electrons have a significant contribution, our data suggest that this stripe antiferromagnetically or¬ 
dered phase in Rbo.sFe 1 . 5 S 2 is a Mott-like insulator with fully localized 3 d electrons and a high-spin 
ground state configuration. Nevertheless, the anisotropic exchange couplings appear to be universal 
in the stripe phase of Fe pnictides and chalcogenides. 

PACS numbers: 25.40.Fq, 75.30.Ds, 75.50.Ee, 78.70.Nx 


Superconductivity emerges in the vicinity of an¬ 
tiferromagnetism (AFM) in both copper based and 
iron based high-transition temperature (high-T c ) 
superconductors [THU] • The AFM in these systems 
share several similarities: antiferromagnetic order in 
layered parent compounds, a spin resonance mode in 
the superconducting state, and the presence of spin 
fluctuations throughout the doping-temperature phase 
diagrams[4]. However, the AFM in the cuprate high-T c 
and iron-based superconductors could have different 
origins. The parent compound of the copper oxide 
superconductors is a Mott insulator with S = 1/2 local 
moments [5]. In the iron pnictides the parent compounds 
are bad metals with several bands crossing the Fermi 
energy. The stripe AF ordering wavevectors coincide 
with the wave vectors connecting the centers of the 
electron and hole Fermi surfaces [6]. In fact many view 
the AF order as due to the Fermi surface nesting. 

From a localized point of view, with 6 electrons in 
the iron 3d orbitals of Fe 2+ , the maximum total spin is 
S = 2. This spin state can be realized when the Hund’s 
Rule coupling energy, Jh , dominates over the crystal- 
field splitting associated with the FeM 4 (M = pnictigens 
or chalcogens) structural unit. On the other hand, a crys¬ 
tal field splitting Acf comparable to the Hund’s coupling 
Jh can lead to an intermediate-spin S = 1 state. In the 
large crystal field extreme, the 3d 6 ions of Fe 2+ will form 
a low-spin singlet S = 0 state[7 9,. In the presence of 


itinerant carriers the spin must be less than S = 2 due 
to charge fluctuations. Thus, while the observation of an 
intermediate-spin state S = 1 does not rule out the pres¬ 
ence of itinerant carriers, the observation of S = 2 would 
require the system to be predominately localized. Not 
surprisingly, the various values of ordered moments ob¬ 
served in different iron-based materials have been inter¬ 
preted in terms of both the local moment picture and the 
itinerant carrier pictnre[7hlQ]. As to the value of the fluc¬ 
tuating local moment, inelastic neutron scattering exper¬ 
iments combined with the moment sum rule revealed an 
increase of S from S « 1 at 10 K to S « 3/2 at 300 K for 
Fei.iTe and a constant S = 1/2 for BaFeoAso pTUT3] . In 
addition, an X-ray emission spectroscopy study was in¬ 
terpreted to imply that the iron spin-state varied between 
S = 0 and 2 in the rare-earth doped Cai_ ;r RE ;E Fe 2 As 2 as 
a function of temperature [14] . In contrast, studies of 
the spin wave excitations in XFe 2 As 2 (X = Ca, Sr, and 
Ba) have found significant special weight contributions 
from both itinerant carriers and local moments [T 3 l fTTR 
EEZ]. These findings suggest that the magnetism of the 
iron pnictides and chalcogenides should be understood 
from a point of view where both itinerant carriers and 
local moments coexist. 

The substitution of sulfur for selenium progressively 
suppresses the superconductivity in Ko. 8 Fe y Se 2 - 2 S^ and 
decreasing the Fe content results in an insulating ground 
state |18j 19 . We found a stripe AF order in insulat- 
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FIG. 1: (color online). Three- (a) and Two- (b) dimensional 
structures of the stripe AF order with rhombic iron vacancy 
order in Rbo.sFe 1 . 5 S 2 . We use the orthorhombic unit cell as 
shown by the solid square in (b) with lattice parameters of 
a ts 5.58 A, b = 5.39 A, and c — 13.889 A. The wave vector 
Q is defined as Q — [H,K,L\ = ( 2nH/a,2nK/b,2nL/c ) in 
reciprocal lattice units (r.l.u). The dashed rectangle is the 
real magnetic unit cell, (c) Dynamic susceptibility x'i^) as 
a function of energy with Ei — 35, 170, and 250 meV at 
8 K. The solid line is computed by the model discussed in 
the text. The dashed line is the dynamic susceptibility of 
BaFe 2 As 2 from Ref. 12 . (d) A candidate for the high-spin 
ground state configuration of the stripe AF order 7. 


ing Rbo. 8 Fe 1 . 5 S 2 with a rhombic iron vacancy order and 
a strikingly similar Neel temperature of T/v = 275 K 
with a moment size of M = 2.8 =b 0.5 /ib as that in 
Kn «iFei ^«Seo [20l |2T] , Photoemission measurements re¬ 
vealed a 980 meV charge gap in the rhombic iron vacancy 
ordered phase [22]. The stripe AF order was proposed as 
a candidate parent compound for the superconducting 
phase in Ao. 8 Fe y Se 2 (A = alkali metal) [20l [23] . It is 
therefore important to characterize the spin waves asso¬ 
ciated with the stripe AF order in insulating Rbo.sFe 1 . 5 S 2 
in order to understand the nature of its magnetism. 

In this paper, we report inelastic neutron scattering 
studies of the spin wave excitations of the stripe AF or¬ 
der in insulating Rbo.sFe 1 . 5 S 2 . Only the spin excitations 
associated with the stripe AF order are observed in our 
experiment, suggesting a nearly 100% stripe AF order 
volume fraction. In the presence of iron vacancy order, 
there are six iron atoms per magnetic unit cell. Hence, 
one expects three doubly-degenerate spin wave branches. 
The first acoustic and the second optical branches are ob¬ 
served clearly in both momentum and energy scans in our 
experiment which uses Ei < 250 meV. The third branch 
is flat in momentum space and can only be observed by 
scans in energy. By fitting the spin excitation spectrum 


to a Heisenberg Hamiltonian with spatially anisotropic 
exchange couplings (SJ\ a = 42 d= 5, SJ\ b = —20 =b 2, 
SJ 2 = 17 ± 2, SJ C = 0.29 =b 0.05 and SJ s = 0.09 =b 0.02 
meV), all of the branches of the spin excitations can be 
accurately described. Furthermore, the total dynamic 
spin fluctuation moment spectrum is calculated to be 
(m) 2 « 20/i^/Fe, similar to that in the block insulating 
AF Rbo.89Fei.5sSe2[24]. After including the contribution 
of the ordered moment 2.8/ie, we estimate the spin to be 
S = 2. Knowing the stripe AF order is an insulator with 
a large charge gap (~ 1 eV), the spin S = 2 suggests that 
all Fe 3d electrons are fully localized. 

Our experiments were carried out on the MAPS time- 
of-flight (TOF) chopper spectrometer at the Rutherford- 
Appleton Laboratory, Didcot, UK, and the BT-7 thermal 
triple-axis spectrometer at the NIST Center for Neutron 
Research, Gaithersburg, USA. We coaligned 1.5 grams 
of single crystals with a mosaic of 1.5° full width at half 
maximum for the two experiments. The stripe AF order 
with M = 2.8 =b 0.5 hb moments and rhombic iron va¬ 
cancy order has been reported elsewhere [21]. For the 
TOF experiment at MAPS, we aligned the c axis of 
the sample parallel to the incident beam at energies of 
Ei = 35, 80,170 and 250 meV at 8 K. The intensities were 
normalized to absolute units by vanadium incoherent 
scattering. For the low energy neutron scattering mea¬ 
surements performed at BT-7, we fixed the final energy 
at 14.7 meV, with horizontal collimations of open-80 / -S'- 
80'-120', where S = sample, and two pyrolytic graphite 
filters after the sample [25.. 

We show spin excitations in the [iL, K] plane at var¬ 
ious energies in Figj2] (a-e). The spin excitations stem 
from the AF wave vectors, disperse outwards and sepa¬ 
rate into two arcs at E = 65 db 5 and 75 =b 5 meV. At 
the energy of 110 ± 9 meV, the wave vectors rotate 90°. 
Weak spin excitations at Q = (±1 d= 0.5,0), (0, ±1 =b 0.5) 
and (±1, ±1) in Fig. [ 2 ] (a) can also be observed. 

To describe the spin waves in Rbo.sFe 1 . 5 S 2 , we em¬ 
ployed a Heisenberg model with in-plane nearest- (Ji a , 
Jib), and next-nearest- (J 2 ) neighbor exchange couplings, 
together with the coupling between layers, J c , as shown 
in Fig. |T] (a) and Fig. 0 (b), and the single ion anisotropy 
term, J s . The Hamiltonian can be written as: 

fl = F^ Sr . s r ,-jT(s“) 2 , (i) 

rr ' r 

where J r y are the effective exchange couplings and (r, r f ) 
label the iron sites [26]. The spin wave excitation spec¬ 
trum can be expressed analytically by solving Eq. (1) 
using the linear spin wave approximation [H3l fT 5 j ITT! l24j . 
We fit the data and convolute the instrumental resolu¬ 
tion using the Tobyfit program[27], From the best fit to 
the experimental data, we determine the parameters as 
SJia = 42 =b 5 , SJ lb = -20 ± 2, and SJ 2 = 17 ±2 meV, 
and for computational convenience an energy indepen¬ 
dent damping T = 7 ± 2 meV. The widths of the spin 
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FIG. 2: (color online). Constant energy slices in the [ H , K] plane of the spin waves at energies of (a) E — 23 ± 3 meV with 
Ei = 80 meV, and (b) E = 45 zb 5, (c) E = 65 d= 5, (d) E — 75 =b 5, (e) E — 110 =b 9 meV with Ei — 250 meV, all at 8 
K. (f-j) Simulations of spin excitations at the identical energies as in (a-e) using the exchange couplings from the best fits to 
the experimental data. The simulations were convoluted with the instrumental resolution. The color bar is the same for each 
energy transfer in units of mbarSr~ 1 meV~ 1 f.uE 1 . 


wave peaks in Ei and K were close to being instrumen¬ 
tal resolution limited as expected for an insulator; this 
also holds true for Ko.8iFei.58Se2[2El GHj- The simula¬ 
tions with the fit parameters at the identical energies of 
Fig. | (-) are presented in Fig. [2] (f-j). 

To compare quantitatively the experimental data with 
the model, we plot cuts along the [if, 0 ] and [1 , if] direc¬ 
tions for a wide range of energies in Fig. [3| The solid 
lines are the results of the best fits with the parameters 
discussed above. The fits are in good agreement with the 
experimental data at all energies. The small discrepancy 
near Q = ( 2 , 0) is due to an acoustic phonon. The weaker 
and flatter cut along the [1 , if] direction at E = 65 =b 5 
meV in Fig. [3] (g) and the cut along the [if, 0] direction 
at E = 110 ±9 meV in Fig. [3] (d) are consistent with the 
splitting of the first branch along the [if, 0 ] direction and 
the 90° rotation of the second branch. 

Fig. [4] (a) and Fig. [4] (b) show the dispersion rela¬ 
tions along the [if, 0] and [1, if] directions with Ei = 250 
meV at 8 K, respectively. The spin excitations from the 
second twin at Q = (0.5, 0),E = 25 meV in Fig. [4] (a) 
and the second branch of spin excitations at energies be¬ 
tween 90 and 120 meV in Fig. [ 4 ] (b) can be observed. 
The dispersion of the spin excitations extracted from ex¬ 
tensive constant energy cuts and Q cuts, together with 
the results of simulations with the best fit parameters are 
plotted in Fig. [4] (c) and Fig. [4](d). Three branches of 
spin excitations can be seen. We tried to fit the disper¬ 
sions in Rbo. 8 Fe 1 . 5 S 2 with the parameters obtained for 
Ko. 8 iFei. 5 sSe 2 ^ j . The dispersions of the first branch 
along the [if, 0 ] and [ 1 , if] directions were matched very 


well, but the second branch along the [l,if] direction 
deviated from the experimental data [see supplementary 
information]. 

In order to determine the exchange coupling between 
layers, J c , and the single-ion anisotropy term, J s , we 
measured the f-modulation of the low energy spin ex¬ 
citations at 2 K. The measurements show that a gap in 
the spin excitations opens up below A = 6 meV and 
that J c only affects the spin excitation spectrum below 
15 meV [Fig. [4] (e)]. By fitting the f-modulated spin 
excitation spectrum, we determined SJ C = 0.29 =b 0.05 
and SJ s = 0.09 =b 0.02 meV. The temperature depen¬ 
dence of the spin gap was also studied and is presented 
in Fig. [4] (f). The spin gap remained sharp right up to 
the phase transition. The scaled magnetic order parame¬ 
ter is plotted along with the temperature dependent spin 
gap. The evolution of the spin gap with temperature fol¬ 
lows the trend of the AF order, in agreement with the 
behavior observed in K 2 N 1 F 4 , a quasi-two-dimensional 
( 2 D) Heisenberg AF insulator [29]. 

To unveil the spin state in the stripe AF order of 
Rbo. 8 Fe 1 . 5 S 2 , we examined the sum rule of the magnetic 
neutron scattering. One can calculate the total fluctuat¬ 
ing moment squared (m 2 ) by integrating the susceptibil¬ 
ity x"(q, lj ) over the band width of the spin excitations 
via 


<m 2 > 


3h r + °° /x /, (q^)dq//dq ^ 
7T Too 1 - exp(—h fjj/k B T ) 


( 2 ) 


The total moment sum rule is Mq = g 2 M 2 + 
(m 2 ) = g 2 S(S + 1 ), where g is the Lande ^-factor and 
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FIG. 3: Constant energy cuts through Q — (1,0) along the 
[. H , 0] (a-d) and [1, K] (e-h) directions at energies of E — 23±3 
meV with Ei = 80 meV, and 45 =b 5, 65 =b 5, and 110 d= 9 meV 
with Ei sss 250 meV, at 8 K. The error bars indicate one 
sigma. The solid lines are the best fits obtained from the 
Tobyfit program. 


M is the static moment. Thus the spin S can be 
extracted [24. 30, 31j. 

The averaged dynamic susceptibility in a Brillouin 
zone x"^) = f x"^ ^)dq/ / dq is plotted in Fig. |l] 
(c). The spin fluctuations in Rbo.sFe1.5S2 are obviously 
stronger than those in BaFe2As2- Integrating the dy¬ 
namic susceptibility through all the spin excitation band 
width results in 29.7±5.5/i^/formula unit (f.u.), and thus 
19.8±3.7/i^/Fe. Taking the ordered moment M = 2.8 ± 
0.5 fis i^o account Mq = (2.8 =b 0.5) 2 ± (19.8 ± 3.7)/i^, 
the total moment squared per Fe is 27.6 ± 4.2/i^, which 
assuming g = 2.0 results in a spin S = 2.2 d= 0.2, which is 
equal to the upper limit of 24/i^ and S = 2 as the Hund’s 
rule result for the 3d Fe 2+ within the error. The results 
reveal that to within the errors all six 3d electrons of Fe 2+ 
are associated with the local moments and in the high- 
spin state. A candidate spin configuration is illustrated 
in Fig. [T] (d). The fact that the carriers are fully local¬ 
ized in Rbo.8Fe1.5S2 is consistent with our photoemission 
measurements on several pieces of single crystals from the 
same batch of nearly 100% stripe AF phase. These mea¬ 
surements also reveal a large charge gap below the Fermi 
energy, suggesting that the stripe AF phase is a Mott- 
like insulator with the integer spin S = 2 23. [32], rather 
than a small gap band insulator [201 12T1 f33j . The rhombic 
iron vacancy order stabilized at a longer in-plane Fe-Fe 
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FIG. 4: (color online), (a) Spin excitations along the [ H , 0] di¬ 
rection, averaging over K = ±0.2 r.I.u and (b) along the [1, K] 
direction, averaging over H = 1±0.2 r.I.u with Ei — 250 meV 
at 8 K. (c, d) The dispersion extracted from experimental data 
and simulations with the best fit parameters. The red circles 
are from the first twin, and the green circles are from the 
second twin [see supplementary information]. The intensity of 
the simulations is proportional to x"(q, a;) x y/E. (e) The 
L-modulation of the low energy spin excitations at 2 K and 
simulations with SJ C = 0.29, SJ s = 0.09 meV and the inten¬ 
sity proportional to x ;, (q, a;). (f) The temperature evolution 
of the spin gap measured at Q = (1,0,1). The solid line is 
the result of a fit to the magnetic order parameter (the blue 
squares) with A( 1 — T/Tn ) 13 , where A is a scaler, T/v = 265 
K, and f3 %= 0.30. 


distance (2.750 A) of the stripe AF order larger than 
that of the block AF order (2.663 A) in Rbo.8Fe1.5S2 
could enhance the correlation and thus promotes the 
localization [23, 132] . 

Several theoretical methods have been successfully ex¬ 
plored to describe the spin waves of the stripe AF order: 
a combination of density functional theory (DFT) and 
dynamic mean field theory (DMFT)f34, 35]; a Heisen¬ 
berg model with the anisotropic in-plane exchange cou¬ 
plings Ji a (> 0), Jib{< 0), and Jo[T3l [T5l [171 [28]: and 
a Heisenberg model with Ji, J2 and a large biquadratic 
coupling Km EH EH- The origins of the anisotropy 
in the Ji a — J15 — J2 model and the large biquadratic 
coupling in the J\ — J 2 — K model are still under 
debate. Long range nematic order is clearly not re- 
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TABLE I: The magnetic exchange couplings and spin states 
in the stripe AF order of iron pnictides and chalcogenides [T3, 

El Eg. 


Compounds 

SJ la 

SJ lb 

SJ 2 (meV) 

S 

M(hb) 

CaFe2As2 

50 zb 10 

-6 ±5 

19 ±4 

1/2 

0.80 

BaFe2 As2 

59 ±2 

-9 ±2 

14 ± 1 

1/2 

0.87 

SrFe 2 As 2 (L) 

31 zb 1 

-5 ±5 

22 zb 1 

0.30 

0.94 

SrFe 2 As 2 (H) 

39 zb 2 

-5 zb 5 

27 zb 1 

0.69 

0.94 

Ko.85Fei.54Se2 

38 =b 7 

— 11 zb 5 

19 zb 2 

- 

2.8 

Rbo. 8 Fe 1 . 5 S 2 

42 zb 5 

-20 zb 2 

17 zb 2 

2 

2.8 zb 0.5 


quired in YFe2As2, since the spin excitation spectrum 
at least at higher energies is little changed at tempera¬ 
tures well above the tetragonal-orthorhombic structure 
transition [T3l H71 ^|39]. The spin waves of Rbo.8Fe1.5S2 
could be described by either model. In particular the 
rhombic iron vacancy order which has already broken 
the C4 symmetry forms at a temperature higher than 
718 K|21]. The anisotropic Ji a and J15 in Rbo.8Fe1.5S2 
could originate from the structural orthorhombicity and 
the possible orbital ordering [40]. For the J\ — J 2 — K 
model, the exchange couplings are estimated to be J\S = 
(Ji a + Ji b )S/2 = 11 =b 3, J 2 S = 17 ± 2, and KS = 
(Jia — Jib)S/A = 15.5 =b 1.4 meV[36]. The biquadratic 
term could be enhanced by the dynamic fluctuations in 
the chalcogen height. Distinguishing the two models mi¬ 
croscopically is beyond the scope of this work. 

We list in Table [I] the fitted magnetic exchange cou¬ 
plings and measured Fe spin values in a number of stripe 
phase Fe arsenides and chalcogenides. The Fermi surfaces 
in these materials vary significantly as do, concomitantly, 
the conductivity, the ordered moments and the effective 
spin values. In spite of this, the exchange couplings mea¬ 
sured in units of SJ are remarkably universal. This result 
is both striking and mysterious. It remains to be seen 
how this relates to the superconductivity in the doped 
materials. 

In summary, we have studied the spin waves of the 
pure stripe AF order in Rbo.sFe1.5S2 over a wide range 
in reciprocal space and energy. Our inelastic neutron 
scattering data reveal that even though the stripe AF 
order has strikingly similar SJ with all the other iron 
pinctides and chalcogenides, it is almost an ideal S = 2 
Heisenberg antiferromagnet with fully localized moments 
inducing Mott insulator behavior. 
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Supplementary: Spin Waves and Spatially Anisotropic 
Electron Interactions in the S = 2 Stripe Antiferromag- 
net Rbo.8Fe1.5S2 
Meng Wang et al. 
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FIG. SI: (color online), (a) Schematic of the iron layer with 
the rhombic iron vacancy order of the first twin, and (b) a 
simulation of the corresponding spin waves at E = 30 meV in 
the [if, K] plane and (c) a simulation of the dispersion along 
the [if, 0] direction with the best fit parameters. The damping 
has been fixed at 3 meV for good viewability. The red arrow 
in (b) indicates the Q in (c). (d) Schematic of the iron layer 
of the second twin, (e, f) The same plots with (b) and (c) for 
the second twin, respectively. 



FIG. S2: (color online). The dispersion extracted from the 
experimental data of Rbo.sFe 1 . 5 S 2 and the simulations with 
the parameters obtained from Ko.siFei. 5 sSe 2 [ 28 ] along the (a) 
[if, 0] and (b) [1 , K] directions. The intensity of the simula¬ 
tions is proportional to x"(q, a;) x VE. 










